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The accuracy of the post-Newtonian waveforms, both in "standard" and "Pade" form, is determined 
by computing their matched- filtering overlap integral with a reference waveform obtained from black- 
hole perturbation theory. 
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I. INTRODUCTION 

The need for very accurate model waveforms to search 
and measure gravitational-wave signals from inspiraling 
compact binaries was recognized less than six years ago 
j|. During this time, post-Newtonian (PN) calculations 
of binary waveforms have been pushed to high order ||^| . 
Currently, the waveform is known to 2.5PN order beyond 
the quadrupole- formula result, and the 3PN waveform 
appears to be within reach. An ever important ques- 
tion is: at which order in the post-Newtonian expansion 
can the waveform be considered to be sufficiently accu- 
rate? In the absence of an exact representation of the 
signal (which might eventually be provided by numerical 
relativity), a definitive answer to this question remains 
elusive. 

An alternative approach to waveform calculations, 
based of black-hole perturbation theory, was also pushed 
forward in the last several years While post- 

Newtonian theory is based on an assumption of slow or- 
bital motion but allows for arbitrary mass ratios, pertur- 
bation theory is based on an assumption of small mass 
ratios but allows for arbitrary velocities. The physical 
situation is that of a small mass orbiting a massive black 
hole, and the emission of gravitational waves is calculated 
by integrating a linear wave equation — the Teukolsky 
equation — in the black-hole spacetime. In the limit 
of a vanishing mass ratio, this approach gives an exact 
representation of the gravitational-wave signal. When 
combined with a slow-motion approximation, perturba- 
tion theory returns an analytic expression for the waves, 
in terms of a series in powers of the orbital velicity. Cur- 
rently, the perturbation-theory waveform is known ex- 
actly in numerical form, and at 5.5PN order in analytic 
form 0. 

Although the perturbation-theory waveform is of re- 
stricted validity, it is still our best guide in trying to 
determine the accuracy of the post-Newtonian expan- 
sion. In this contribution to these proceedings, this issue 
will be examined in two different contexts. First, sup- 
posing that gravitational-wave signals from small-mass- 
ratio binaries can be picked up by such interferometric 



detectors as LIGO/VIRGO, I calculate (Sec. 3) the loss 
in signal-to-noise ratio incurred by using second post- 
Newtonian waveforms as search templates. My results 
suggest that the 2PN waveforms will be acceptably ef- 
fective for searches. Second, I review (Sec. 4) the very 
interesting proposal by Damour, Iyer, and Sathyaprakash 
|^| to use Pade approximants to improve the convergence 
of the post-Newtonian series. My results suggest that 
a Pade version of the 3PN waveforms (yet to be pro- 
duced) will make very accurate templates for searches, 
and possibly also for measurements. In Sec. 2, I present 
the framework used throughout this contribution. 

The work presented here was partially carried out with 
Serge Droz. Additional details can be found in MM. 



II. PRELIMINARIES 

The output of an idealized gravitational-wave detec- 
tor is written as o(t) = n(t) + s(t), where n(t) is sta- 
tionary Gaussian noise and s(t) is the signal. The de- 
tector output is analyzed by matched filtering against 
a bank of templates h(t;X, 0), Here, A collectively de- 
notes the template's kinematical parameters (arrival time 
and initial phase), while 8 collectively represents the 
dynamical parameters (the two masses, assuming that 
the binary companions are nonrotating) . The statisti- 
cal properties of the detector noise are summarized by 
its spectral density, denoted S n (f), where / is the fre- 
quency. Throughout we use the useful approximation 0] 
S n (f) = 5 e(/-/ min )[(/ //) 4 + 2 + 2(/// ) 2 ], where the 
parameters So, /mm, and /o are set by the detector. For 
the initial LIGO, / min = 40 Hz and f = 200 Hz, while 
./mm = 10 Hz and fo = 70 Hz for the advanced LIGO. 
The value of Sq is irrelevant for our purposes. 

We define the ambiguity function A(\, 6) by 



A(X,9) = 



(s\h) 



where, for any functions a(t), b(t), 



(a\b) = 2 



a*(f)kf) + a(f)b*(f) 

Sn(f) 



df 



(1) 



(2) 



1 



is an inner product in the Hilbert space of gravitational- 
wave signals (H). Here, a(f) is the Fourier transform of 
a(t), and an asterisk denotes complex conjugation. The 
ambiguity function measures the Hilbert-space "angle" 
between the signal and the template. This "angle" varies 
with the template parameters, and we are interested in 
maximizing the ambibuity function over these parame- 
ters. We construct the semi-maximized ambiguity func- 
tion (SMAF) by maximizing only over the kinematical 
parameters: 



SMAF(6>) = m&xA(X,e). 



(3) 



The SMAF will be computed in Sec. 4. The fully max- 
imized ambiguity function is Apostolatos' fitting factor 



FF = max.4(A,0). 

. A 



(4) 



The fitting factor is equal to ratio of the actual signal-to- 
noise ratio, obtained with an imperfect set of templates, 
to the signal-to-noise ratio that would be obtained if a 
perfect set of templates were available. Consequently, 
the loss in event rate due to template imperfection is 
1 — FF 3 . The fitting factor will be computed in Sec. 3. 

In the calculations of Sees. 3 and 4, the reference sig- 
nal will be provided by black-hole perturbation theory, as 
was discussed in Sec. 1. Keeping only the dominant mode 
at twice the orbital frequency, the frequency-domain sig- 
nal can be expressed as M 



5(f) a: A(v) e l * M 



(5) 



where v = (nMf) 1 / 3 is the orbital velocity (M is the 
total mass of the binary system, and /i will denote the 
reduced mass); A(v) is the signal's amplitude, and ^(v) 
its phase. When using the stationary-phase approxima- 
tion to calculate the Fourier transform, the phase is found 
to be given by 



T , . 5M 



" (v 3 - v' 3 )Q(v') 
v /9 P(v') 



dv'. 



(6) 



Here, P(v) = (dE / dt) / (dE / dt)op is the rate at which 
the gravitational waves remove orbital energy from the 
system, normalized to the quadrupole-formula expres- 
sion (dE/dt) QF = -(32/5)(/VM)V°; and Q(v) = 
(dE/dv)/{dE/dv) N = (1 - 6v 2 )/{l - 3v 2 ) 3 / 2 gives 
the differential relation between orbital energy and or- 
bital velocity, normalized to the Newtonian expression 
(dE/dv)-^ = —/J.V. The functions A(v) and P(v) are cal- 
culated by numerically integrating the Teukolsky equa- 
tion for a small mass moving on a circular orbit of the 
Schwarzschild black hole. 

On the other hand, the templates are given by 



h(f;\,6)<xv 



-7/2 e ii>{v;X,0) 



(7) 



where v 1 1 2 is the leading-order approximation to A(v), 
and ip(v; A, 6) is a post-Newtonian approximation to 



*(«). In Sec. 3, ip(v;A,0) will be given by the 2PN 
approximation obtained by Blanchet, Iyer, Will, and 
Wiseman M using post-Newtonian theory. In Sec. 4, 
?/>(i>; A, 0) will be given by the slow-motion, perturbation- 
theory results of Tanaka, Tagoshi, and Sasaki JlTJ]; this 
approximation does not incorporate the finite-mass-ratio 
terms of the true post-Newtonian series, but it is carried 
out to a high PN order. 

It should be emphasized that the perturbation-theory 
results are valid only in the limit fi/M — ► 0. Neverthe- 
less, I will still consider binary systems with large mass 
ratios, so as to give an indication of the robustness of my 
conclusions. 



III. SECOND POST-NEWTONIAN WAVEFORMS 
AS SEARCH TEMPLATES 

Serge Droz and I have calculated fitting factors for sev- 
eral binary systems using 2PN waveforms as templates. 
Our results are presented in the table, where for com- 
parison I also display the fitting factors for Newtonian 
templates. 





advanced LIGO 


initial LIGO 


System 


FF-Newton FF-2PN 


FF-2PN 


1.4+1.4 M Q 


79.2% 


93.0% 


95.8% 


0.5 + 5.0 M Q 


51.6% 


95.2% 


89.7% 


1.4+10 M Q 


55.8% 


91.4% 


88.4% 


10 + 10 M 


70.1% 


91.7% 


86.3% 


4 + 30 M Q 


61.3% 


86.6% 


67.9% 



We see that for the advanced version of LIGO and for 
most of the binary systems considered, the fitting factors 
are all above the 90% mark. According to Apostolatos' 
criterion second post-Newtonian waveforms should 
make acceptably effective search templates. Additional 
details can be found in Ref. M . 



IV. PADE APPROXIMANTS 

Recently, Damour, Iyer, and Sathyaprakash (DIS) || 
discovered a clever way of expressing the post-Newtonian 
waveforms in terms of Pade approximants, which acceler- 
ates the convergence of the post-Newtonian series. Here 
I shall restrict all considerations to the test-mass limit, 
fj,/M -> 0. 

It is known that Pade approximants work especially 
well for rational functions. However, the orbital energy 
(per unit rest mass) of a test particle in Schwarzschild 
spacetime is not a rational function of the orbital ve- 
locity: E(v) = (1 - 2u 2 )/(l - 3u 2 ) 1/2 . To circumvent 
this, DIS choose to employ the alternative energy func- 
tion e(v) = E 2 - 1 = -u 2 (l - 4u 2 )/(l - 3w 2 ), which is 
rational. In a post-Newtonian context, only the first few 
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terms of a Taylor series about v = can be calculated. 
For example, at the 2PN level, e(v) ~ -v 2 (l - v 2 + 3v 4 ). 
It is easy to check that the equivalent Pade approximant 
coincides with the exact result. 

The Pade series therefore converges rapidly to the ex- 
act result for e(v), and it may be used to estimate the 
position of the pole: v po \ c = In the test-mass 

limit, this estimate is exact. It is possible, however, that 
the estimate might not be as accurate in the finite-mass- 
ratio case, for which an exact expression is not available. 
As a way of challenging the DIS proposal, I will allow 
for a possible uncertainty in the position of the pole, by 
introducing a free parameter £ such that v po i e = £/V3- 

The second element of the DIS proposal is to express 
P(v), the luminosity function introduced in Sec. 2, also 
in terms of Pade approximants. This requires some 
thought, because schematically, the post-Newtonian se- 
ries for P(v) (as calculated by Tanaka, Tagoshi, and 
Sasaki (Tol) takes the form 

P{v) = 1 + v 2 + v 3 + v 4 + v 5 + (1 + lnv)v 6 
+ v 7 + (1 + lnv)v s + (1 + lni>)i; 9 
+ (1 + bxv)v 10 + (1 + \nv)v 11 ; (8) 

the presence of logarithmic terms prevents a straightfor- 
ward conversion to a Pade form. The way around this, 
as implemented by DIS, is to factorize these terms. DIS 
also choose to factorize a simple pole at v = v po \ c , be- 
cause this pole, already present in the energy function, 
can be shown to appear also in the luminosity function. 
The new expression for P(v) is therefore 

P(v) = [l + {v e + --- + v n )lnv] 

X [l - v/v po i c ] 1 

x [l + v + ••• + «"], (9) 

and DIS re-express the third factor in terms of its equiv- 
alent Pade approximant. 

Using the DIS proposal, I have calculated semi- 
maximized ambiguity functions (SMAF) for selected bi- 
nary systems. The results depend on the order n at which 
the post-Newtonian expansion is truncated. [For exam- 
ple, n = 6 corresponds to Eq. (^) with all terms of order 
v 7 and higher discarded; this is the 3PN-Pade approxima- 
tion.] The results depend also on £, which parameterizes 
the possible uncertainty in the position of the pole. The 
values {1.0,0.9, 1.1, oo} were selected; the choice £ = oo 
corresponds to the absence of the second factor in Eq. (|]) 
— the pole is not factorized. 

The representative case of a 1.4+10 M© binary system 
is represented in the figure, which plots the SMAF as a 
function of the truncation order n, calculated using a 
noise curve appropriate for an advanced LIGO detector. 
The curve labelled "T" is obtained by using the usual 
Taylor expansion (|^) for the luminosity function, together 
with a Taylor expansion for Q(v). The curves labelled by 
numbers are obtained by using the Pade form of Eq. (||), 




k' ' 5 ' ' t 7 ' ' ' ' h b' ' ' 10 ' Tl 

n 



with £ given by the corresponding number. The curve 
labelled "oo" is also obtained by using the Pade form of 
Eq. (g), but in the absence of the second factor. The 
exact expression for Q(v) is used for all the Pade curves. 

We see that the Pade curves converge to unity much 
more rapidly than the Taylor curve. This demonstrates 
the great success of the DIS proposal. Furthermore, this 
is true for all the Pade curves, whether or not they in- 
corporate an uncertainty in v po i c . I find this robustness 
of the DIS method quite remarkable. It is also interest- 
ing to note that the choice £ = 0.9 produces the largest 
values for the SMAF. 

These results suggest that 3PN-Pade waveforms will 
make a suitable choice of search templates. They may 
also be sufficiently accurate for the reliable estimation of 
source parameters. Additional details will be presented 
elsewhere ||. 
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